Quantum heuristic algorithm for traveling salesman problem 
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We propose a quantum heuristic algorithm to solve a traveling salesman problem by generalizing 
Grover search. Sufficient conditions are derived to greatly enhance the probability of finding the 
tours with the cheapest costs reaching almost to unity. These conditions are characterized by 
statistical properties of tour costs and shown to be automatically satisfied in the large number limit 
of cities. In particular for a continuous distribution of the tours along the cost we show that the 
quantum heuristic algorithm exhibits the quadratic speedup over its classical heuristic algorithm. 



I. INTRODUCTION 

The quantum computation are potentially more power- 
ful than its classical counterparts [lj. In particular some 
quantum algorithms enabling the dramatic speedup has 
opened new possiblities of quantum computation to solve 
hard problems R is usually considered that the 

hard problems have so many candidate solutions that it 
is hard (or practically impossible) to search the exact 
solution among them quickly. Such hard problems are 
usually classified as "non-deterministic polynomial (NP) 
class" in complexity theory @ . Thus researchers in this 
field have expected a succession of quantum agorithm to 
solve the NP problems, depending on quantum superpo- 
sition which enables to investigate simultaneously all the 
states. However the development of quantum algorithms 
has been rather stagnant, with any remarkable quantum 
algorithms having not been discovered in the last few 
years 3 . and it is still unclear that all NP problems can 
efficiently be solved by quantum algorithm. 

On the other hand, in computational science, heuristic 
approach is one of the most efficient ways to solve the NP 
problems. In heuristic approach, the problem is usually 
simplified at the cost of accuracy, and thus it may fail to 
find the exact solution [8]. Heuristic algorithm, instead, 
enable us to gain computational speedeup by finding the 
approximate solution, near to the exact one, so quickly. 
Therefore it is natural that current research in quantum 
algorithm has been directed towards using this heuristic 
approach However, it also has not been elucidated 

thoroughly that the heuristic approach is always helpful 
for speedup when we deal with the problem in a quantum- 
mechanical way. The present work is actually motivated 
by the goal of investigating this issue. 

As aforementioned, it is generally known that the 
heuristic approach is useful to solve the hard problems. 
However, in the case of using the quantum system, this 
remains valid only for the search problem [9|, [ll|. In 
addition, the heuristic algorithms are problem-specific: 
There is no guarantee that a heuristic procedure used for 
one problem will be applicable to others (l2| - tT3 |. There- 
fore we consider a particular set of problems, called "NP- 
complete" , because of the property that all NP problems 
can efficiently be mapped into a problem in this class Q. 



Typical examples belonging to the NP-complete class 
are a traveling salesman problem (TSP) and an exact 
cover [isj] . We consider exclusively the TSP in this pa- 
per. Thus, if a effective heuristic algorithm for TSP were 
found, it could be mapped into a procedure for solving 
all NP problems even though problem-specific. 

In TSP, a salesman travels n cities by visiting each 
city once and only once, and returning to the starting 
city. The problem is to find the lowest-cost tour among 
(n — 1)! possible tours, where the symbol "!" stands for 
the factorial. Throughout the paper, we assume that the 
oracle answers only the overall cost for the given tour 
when queried, hiding the details on the costs of city pairs. 
This assumption is made in order to compare the per- 
formances of algorithms under the equal conditions, ne- 
glecting any improvement by structural information (see 
Ref. [161 ] for such an improvement). With such an or- 
acle of the tour costs, a classical exact algorithm is to 
query the costs of all the tours until finding the cheap- 
est. As the cheapest cost is unknown, contrary to the 
data search where the target is known Q, the search 
should be done over all the tours. This demands a huge 
number of (n— 1)! ~ 0(e n +n ) queries, independent of 
whether the algorithm is deterministic or probabilistic. 

Many quantum approaches for the TSP (or equivalent 
to this problem) have been studied @, [l(| fl7l - [l9} . In par- 
ticular the quan tum annealing method was investigated 
for the TSP [ldl.[l9j. It has been discussed, however, that 
finding a ground state becomes extremely difficult as in- 
creasing the size of system (or equivalently the number of 
cities) [H, [IH ■ R is unclear that the quantum annealing 
can be applied to the extremely large number of cities. 
On the other hand, one may convert a TSP to the prob- 
lem of finding a minimum for a given cost function, called 
"a minimum search." A quantum algorithm to obtain an 
exact solution of minimum search was proposed, based 
on Grover search, which works in 0(V~N) queries, where 
N is the number of elements (or tours in TSP) 20]. This 
is the quadratic speedup of the classical exact algorithm 
earlier, as \/~N ~ S 2 nlnn +2 n . 

We instead consider a classical heuristic algorithm 
which is adapted to the circumstance: The tours along 
the cost c arc distributed, close to a normal distribution 
with the average c and the standard deviation Ac in the 
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limit of the large number of cities, n —¥ oo. This circum- 
stance has commonly been considered in the literatures 
[1, EH- (The statistical information can easily be esti- 
mated by sampling a small number of tours.) Taking an 
approximate solution to be a tour whose cost is less than 
Ac, the probability / of finding a solution is asymptoti- 
cally given, 



■ In n 



(i) 



Then the classical heuristic algorithm is following: Query 
1// tours that are randomly chosen and find the cheapest 
one among those tours. Then one solution is "likely" to 
be found within 0(1/ f) number of queries. Note that 
this classical heuristic algorithm is much faster than the 
exact quantum algorithm by removing the logarithm in 
the exponent. 

In this paper, we propose a quantum heuristic algo- 
rithm for the TSP by generalizing Grover operation. Suf- 
ficient conditions are derived under which the algorithm 
can greatly enhance the probability of finding approxi- 
mate tours with the cheapest costs, reaching almost to 
unity. These conditions are characterized by statisti- 
cal properties of tour costs and shown to be automat- 
ically satisfied in the large number limit of cities. The 
present quantum heuristic algorithm is shown to work in 
0(1/ '\/J) queries, a quadratic speed up of the classical 
heuristic algorithm with its order 0(1/ f) in Eq. (JTJ. 



II. QUANTUM TRAVELING SALESMAN 
PROBLEM 



To represent the tours quantum-mechanically, an or- 
thonormal basis set of N states are appropriate and they 
are fully discriminated by a von-Neumann orthogonal 
measurement. We assume that a quantum register is of 
N dimension for a sake of simplicity and its basis states 
\Ti) respectively correspond to the tours 

An oracle operation C is introduced to answer the tour 
cost, when applied, such that 



C\T) 



\T), 



(2) 



where <p(T) stands for the overall cost, given tour T, and 
it is called a cost phase. Here, every cost phase is defined 
in between and 2n by scaling the tour costs. The scaling 
can easily be made once obtaining the strictly lower and 
upper bounds of the tour costs [lfj ■ It is remarkable that 
the oracle constructed will not expose any details such as 
the city-pair costs but just the overall tour costs. In order 
to construct the cost oracle, it is convenient to define a 
tour state as \T) = |T) = |C )®|Ci)®- • ■^\C n -x), where 
|Cfe) stands for the k-th visiting city, and \Cq) is both the 
starting city and finally visiting city [2l[ . The conditional 
operation Pjk is also introduced for cost oracle, such that 
Pjk \Cj) \Ck) — e lc i k \Cj) \Ck), where Cjk denotes the cost 
from j-th to k-th visiting city. Since this conditional 
operation Pjk is efficiently be decomposable to two-qubit 
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FIG. 1: Simple quantum circuit for cost oracle. The box 
coupled with another channel stands for the conditional phase 
operation Pjk. The cost oracle can efficiently be constructed 
with n number of conditional phase operations. 



and/or Hadamard gates [22|], we can finally construct a 
cost oracle C with relatively small computational efforts 
by using the Pjk (as depicted in Fig. [1]). 

III. GENERALIZED GROVER OPERATION 

We generalize Grover search Q and apply its general- 
ization to amplifying the probability for approximate so- 
lutions. Before the generalization, we briefly discuss the 
Grover operator Go in an alternative perspective. Apply- 
ing Go iteratively leads an initial state |i/ , o) to a target 
IV'i). [G.l] It is a composition of two reflections: 



Go — I^ I^ 1} 



(3) 



where 7^ \ipi) = - and J^ 4 \ip±) = \ipt), V^) 
orthogonal to \ipi). [G.2] Small overlap between ipf. 
\/J <C 1, where \f] = (^ol^i) is assumed to be a real 
number for a sake of simplicity. [G.3] Go is a rotation 
of angle 2yf] except a global phase, up to the first order 
of y/J, on the plane spanned by |^o) and Letting 

\4>£) = {M-V7\1>i))/VT zr I, 



Go 



■exp(i2y/fa y )®i c 



(4) 



where 1 is an identity operator, a y a Pauli spin oper- 
ator in the subspace spanned by {l^i) ; |V , f L )} I an d l c 
is an identity operator defined in the complementary 
subspace. Note that |^o) — e '*v(iVJ ~y) ^i)- [G.4] 
The number of iterations from (^o) to the target is 
R = [(tt/2 — y/J)/2y/~f \ , where [x] is the nearest integer 

to x, that is, Gq |^o) — l^i)- 

We propose a generalized Grover operator, defined by 



G 



-Iij) C, 



(5) 



where G is the cost oracle in Eq. ([2]) and I^ \i/jq) = — \i/)q) 
while leaving other orthogonal states invariant. Here, the 
initially prepared state \ipo) of the quantum register is 
given by 



(6) 
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where the sum is over the all N of tours. The generalized 
Grover operator G in Eq. (|5|) consists of not two reflec- 
tions but a single reflection and the cost oracle. This is 
the distinct feature of the generalized Grover operator, 
different from those investigated in Ref . [23|, [24| • 

The definition of G is made by observing two simple 
cases. In each case, the tours are divided into several 
groups according to their costs. The j-th group contains 
the number Nj (c) of tours and their costs c are all to be 
encoded to a certain cost phase <f>j. The first case is two 
groups: One of iVi tours having <\>\ — 0, and the other of 
Ni = N — Ni tours having <p2 = 7r. Assuming Ni <C N, 
G is directly approximated to the G rover operator Go, 
upto the first order of y/J = y/Ni/N, We can then find 
the tours of Ni group with near unity probability. 

The other simple case is that the tours are divided into 
four groups, each of which has the cost phase <fij — jir/2 
for j = 0, 1, 2, 3, respectively. Assume further that N ~ 
N 2 ~> Nj for j 7^ 2. It is convenient to define four states 

\^) by 

|^> = ^=V|T>, (7) 



where Sj denotes the j-th group, and hence the sum 
is over the tours T in Sj. The four states form an or- 
thonormal basis {l^j)}, which spans a four dimensional 
subspace. In particular, the initial state \ipo) is expanded 
in terms of the basis as 



It can be seen that G is a transformation between \tt) 
and the other basis states. On the other hand, the fourth 
power 



G 4 fa exp(i8 v /7o(T 2/ ) © i c 



(12) 



is a rotation of angle on the plane spanned by \ir) 

and |0), leaving the other two states \ir/2) and |37r/2) 
unchanged. Here a y is a Pauli spin operator in 1 0) - 1 vr) 
subspace and l c is an identity operator in |7r/2)-|37r/2) 
subspace. It is remarkable that, in our approach, we do 
not have to generalize the "phase matching condition" 
25]. Noting |Vo) « V^k), then (tt/2 - V7o)/8V7o 
iterations of G 4 approximately transform 



— ► y/hThH 



/ Ak/2) + 
VAk/2> 



T*W) + V%|3tt/2) 
Vh\^/2). (13) 



IV'o) = £ y/fj \<t>i) 



(8) 



It is clearly seen that the probability of finding |0) is 
greatly enhanced as / + /2 — 1 ^> A.3, reminding of 
H ^> /o)/i,/3- We note that the number of iterations 
needs to be a multiple of 4 to exclude the possibility of 
leaking the probability to the other tours. One might be 
concerned about the case that /i,3 > /q but this does not 
significantly alter our main result as /i,3 <C f 2 - 

This procedure is generalized straightforwardly to 2M 
groups of tours, j-th group having phase <f>j — jn/M, 
where j = 0, 1, 2M — 1. The arguments similar to the 
case of 4 groups can be applied now with 2M dimensional 
subspace. Then, the 2M-th power of G is given by, up 
to the first order of Wfj, 



where fj — Nj/N. The subspace spanned by {|</>j)} re- 
mains invariant under the transformations by the cost 
oracle C and the reflection I^ : C \<j)j) = e l ^ j \<f>j) and 

ko \<t>j) = l^i> " 2 V N i/ N IV'o) = Ej c 3 Wj) with some 
real numbers Cj. As long as staying in the subspace, 
both operators can be represented in terms of In 
particular, the cost oracle C is reduced and equivalent to 



C s 



-1, 



Ed 



\<t>j) (<f>j 



(9) 



where l s = Ej l&j) ('Pjl ^ s an identity operator in the 
subspace. The necessary condition similar to [G.2] is cast 
in the form of 



(Vol (C s + i.) \iPo) 



0, 



(10) 



which is equivalent to yffj -C 1 for j 7^ 2. We may 

regard C s as the cost oracle, omitting the subscript s, 
without loss of any generality. Within the subspace, the 
generalized Grover operator is approximated, upto the 
first order of y/J] for j ^ 2, as 



G « - |0> <0| 
3 



i|7r/2) (?r/2 1 +i|3vr/2) (3tt/2| 



-2^75 (lW2) (tt| - e^ 2 k) 0V2|) • (11) 



G 



1M 



exp(i4M yffady) © i c 



(14) 



where Efc exp(i2irkj /2M) = 2M5jo is used and & y is 
defined in |0)- |vr) subspace, the same as the cases of M = 
1,2. In order to obtain the almost- unit probability of 
finding |0), the number of iterations of G 2M is given by 
(tt/2 - v / 7o)/4-M"v / 7o : The number of iterations R of G 
or the number of queries becomes 



R = 



2V7o 



(15) 



j¥2 



Here an accumulated error e for very large R is negligible 
as given by e - 0(y%). 

We analyzed that the present algorithm works in which 
the discretization of the tours into 2M groups is admit- 
ted and N(ir) N(<f>j) for all the other phases cf>. In 
other words, if the cost phases <pj = and <f>j = tt are 
well defined and their tour population ratio fo / f n is much 
less than unity, the presented algorithm can be applied to 
solving TSP. 



IV. VALIDITY FOR CONTINUOUS MODEL 

We investigate applicability of our algorithm in a con- 
tinuous tour-cost model with its statistical properties. In 
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other to calculate the total iterations, a Gaussian model 
is used as a typical example @, UU . To begin with, we 
define the density of tours v(c) taking a given cost c, in- 
stead of the number of tours JVj (c) in a discretized j-th 
group. The density of tours v(c) is defined, in the con- 
tinuum limit of n — > oo, by 



i/(c)=£*(c(T) 



(16) 



where c(T) is the overall cost of the given tour T and 
5(c) is a Dirac delta function. We note here that in the 
continuum limit the initial state \ipo) is rewritten as 



lV>o) 



/o v JV 

where the states |^>) is defined by 



* £^-0(T))|T) 



(17) 



(18) 



Such a definition of \<p) is originated from the Eq. (J7J. 
Here, these states |(/>) form the unnormalized basis, sat- 
isfying <<^> = 5(f - <t>). 

Then, we investigate the statistical properties of tour 
costs, which are cost average c and standard deviation 
Ac. The c and Ac are easily estimated, assuming that 
city-pair costs Cjk are randomly chosen in the interval of 
[01,02] with an identical probability density distribution 
P(cjk) — l/(c2 — ci). A tour cost c(T) of n cities is in 
the interval [nc\ , nc 2 ] ■ The c is given by 



1 



(n-1) 



n— 1 n— 1 

E^) = ^EE^ 



(19) 



j=0 fe#j 



Here the last equality holds as every Cjk appears (n — 2)! 
times in all tours and c(T) = (n — 2)! Cjfc, where 
the sum runs over k ^ j. Noting that the ensemble 
average of every single pair cost Cjk is equal to (c 2 — ci)/2, 
the cost average c is estimated to 



c C st = n- 



C2 - Ci 



(20) 



The average of the second moment is given by 



r-2 - 



(n-1)! 



T 

n-1 



(n-2)\ v 2 
(n-1)! I-*' J1J2 



2(n- 3)! 

( n _ 1)1 ^71. 



E 



("-3)! 
(n-1)! 



E 



(21) 



where the last equality holds similarly to Eq. (|T9|) . The 
ensemble average of every c 2 ^ is equal to (c^ — cf )/3(c2 — 
ci). On the other hand, every Cj 1 j 2 Cj 3 j 4 is averaged to the 



simple product of the two ensemble averages of Cj U2 and 
c hji , i.e. [(c 2 - Ci)/2] 2 , (ji,j 2 ) 7^ (J3,J4)- The second 
moment is estimated to 



C2 + C 2 Cl 
C est = 77- 



(ri 2 — n) 



(c 2 - ci) 



3 ' v " "' 4 

The estimated variance is then given as 



Ac 2 f = c 2 



10c 2 ci 



12 



(22) 



(23) 



Then the statistical properties of tour costs, the Eq. (|2"0")) 
and Eq. (|23|) . are briefly represented as 



Ac 

c 



1 



(24) 



The distribution u(c) is transformed to v(4>) with respect 
to the cost phase 4>, by the relation <fi(T) = 2irc(T) / (7ic 2 — 
nci). It is notable that the average cost phase cf> = it 
and the standard deviation Atp ~ 1/ y/n. In the limit of 
71 — y oo, the tours population becomes concentrated to the 
average <fi = n as A<p — > 0. 

We shall finally investigate the applicability of our 
algorithm to the continuous tour-cost distribution of 
n — > oo. Remind that in our approach the crucial condi- 
tions are that (a) the two cost phases of = and (f) = it 
are well defined and (b) their population ratio fo/ f« is 
much less than unity. The condition (a) is rephrased as 
the existence of 77 <C 1 such that 



C 



n/2 



n/2 



v{4>) 

N 



C 



■n/2 



n/2 



N 



TolO) (25) 



e i7r V/^k), (26) 



where C is the cost oracle. The cost phase = 0, de- 
fined in Eq. (|23|) . does not distinguish the lowest costs 
(from to 77/2) and the highest costs (from —77/2 to 0). 
Heuristic algorithms do not in general provide the indi- 
cations to verify that an obtained solution is indeed op- 
timal. Nevertheless we can discriminate the cost phases 
based on quantum phase estimation technique with an 
ancillary register to measure the estimated cost phase 
[Hj]. If A(f> -C 77 < 1, it is clearly seen that Eq. (|26]l 
holds as v{4>) ~ N5(<f> — tt). For the phase 0, noting 
77 <C 1, we expand C\<p) in Eq. ([23)1 upto the first order, 
i.e., C\4>) ~ \4>) + i<fi\(f>). The first order is negligible as 



-n/2 



12 



(27) 



Thus, the condition (a) holds for both phases and 7r if 
A<p By using Eq. (JMJ), the sufficient condition 

is rewritten as 



77 1 
1 <C -7^7 *C -7—7 

A(j> Ac/) 



(28) 
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This is also sufficient for the condition (b), as the most of 
tours are concentrated around the average <f> — tt within a 
few of A(j>. Thus the conditions (a) and (b) are character- 
ized by the statistical properties in a continuous tour-cost 
model. 

To calculate the total iterations, we consider a Gaus- 
sian density of tours v(c) near the lowest costs 



u{c) 



AT V ° 

N—— exp 
Ac 



{c-c? 
2Ac 2 



(29) 



where c is the average of tour costs and Ac the stan- 
dard deviation: c = J dc cv(c)/N and Ac 2 = J dc (c — 
c) 2 v(c)/N . The normalization factor vq is determined 
such that J dc v{c) — N, the total number of tours. Then 
the population fo is given by 



fo- W 



n/2 



77/2 



dip v(4>) ~ e 



(30) 



This is consistent to Eq. ([T]) for the classical heuristic 
algorithm. Nevertheless, the number of queries in our 
algorithm is O(l/v / /o), as seen in Eq. U5\) . leading to 
the quadratic speedup of the classical heuristic algorithm. 



V. SUMMARY 

As the classical heuristic algorithm is much faster than 
the exact quantum algorithm, we proposed a quantum 



heuristic algorithm to solve the TSP by generalizing the 
structure of Grover search. The generalized operation 
consists of the conventional reflection with respect to the 
equal superposition state and the cost oracle which an- 
swers the cost phase when queried, distinct from the pre- 
vious generalizations. Sufficient conditions were derived 
to greatly enhance the probability of finding the cheap- 
est tours and they were characterized by the statistical 
properties of tour costs. In particular for the Gaussian 
distribution of the tours along the cost, we showed that 
the conditions are satisfied and the proposed quantum 
heuristic algorithm exhibits the quadratic speedup of its 
classical heuristic algorithm. It is thus an open ques- 
tion whether the quadratic or faster speedup can still be 
achieved when employing the structural or geometrical 
information of the city-pair costs. 
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